基于NPSVM模型的超新星识别方法

巡天观测与高能物理、黑洞天文等领域均有密切的联系.基于星系-超新星二分类问题,研究光谱数据预处理,结合余弦相似度改善PCA(Principal Component Analysis)光谱分解特征提取方法,用SDSS(the Sloan Digital Sky Survey)、WISeREP(the Weizmann Interactive Supernova data REPository)组成的5620条光谱数据集训练支持向量机,可以得到0.498%泛化误差的识别模型和新样本分类概率.使用Neyman-Pearson决策方法建立NPSVM(Neyman-Pearson Support Vector Machine)模型可进一步降低超新星的漏判率.     

基于双谱滤波的综合脉冲星时算法研究

毫秒脉冲星的自转频率非常稳定,提供了一种独立的基于遥远自然天体并能持续数百万乃至数十亿年的时间基准,具有稳定性强、运行时间长、服务范围广等特点.为了减弱毫秒脉冲星计时观测中各种高斯噪声对脉冲星时的影响,研究了一种基于双谱滤波的综合脉冲星时构建算法,处理分析了国际脉冲星计时阵(International Pul-sar Timing Array,IPTA)最新发布的4颗毫秒脉冲星(PSR J0437-4715、J0613-0200、J1713+0747和J1909-3744)的观测数据,分析了不同时间尺度综合脉冲星时的稳定性,并与构成国际原子时(International Atomic Time,TAI)的4家授时单位原子钟稳定性进行了比较.结果表明:双谱滤波算法能够较好地抑制观测噪声,提高综合脉冲星时的稳定性.相比于经典加权算法,综合脉冲星时1 yr、10 yr稳定度从7.77×10^-14、8.56×10^-16分别提高到1.50×10^-14、3.50×10^-16,单脉冲星时稳定性的提升也类似.同时发现,综合脉冲星时稳定性在5 yr及以上时间尺度上优于原子钟稳定性,可用于改善当前原子时的长期稳定性.     

黄土洞室掘进深度与围岩位移变化规律研究

应用PLAXIS 2D/3D有限元分析软件,对某黄土洞室进行了二维、三维分步掘进、三维一次掘进等工况的弹塑性有限元分析,研究分析了不同工况下洞室围岩位移的变化规律,并对掌子面的空间效应进行了分析。分析结果表明,二维与三维分步计算的最终位移基本一致,可以用二维的计算结果来估算洞室开挖的最终位移。三维分步掘进的分析表明,不同埋深的洞室其位移比的变化规律基本一致,当L/B<1时(L为掘进深度,B为洞室跨度),位移比呈线性增加;当L/B>1时,位移比呈非线性增加;当L/B>3时,其最大位移接近最终位移。掌子面的空间效应分析表明,在掌子面处的位移约为最终位移的1/3,距掌子面0.5倍跨度处的位移约为最终位移的2/3,距掌子面2.5倍跨度处的位移基本达到了最终位移。     

基于GPS和水准数据的微粒群算法联合反演祁连山北缘断裂三维滑动速率

将微粒群算法与位错理论模型相结合，采用中国地壳运动观测网络提供的青藏高原地区2001~2004年GPS测量数据和2000~2006年水准测量数据，通过常规定权和附有相对权比的方法对祁连山北缘断裂的三维滑动速率进行联合反演，并与蚁群算法反演结果进行对比。结果表明，微粒群算法收敛速度快、稳定性高，结合经典位错理论模型，是一种可以有效求解断层三维滑动速率反演问题的优化算法，在大地测量反演领域极具应用潜力。     

微量硫化物、硫酸盐硫同位素EA-IRMS在线测试分析

常规EA-IRMS硫同位素测试中, 硫化银(Ag2S)的需样量为0.2~1.0 mg, 硫酸钡(BaSO4)的需样量为0.35~1.5 mg, 较大的需样量已难以满足珍贵样品及微区样品的分析要求, 因此如何减少测试样品量已成为EA-IRMS测试分析工作中急需解决的问题。通过对EA-IRMS测试系统的研究发现, 在采用常规方法测试样品时, 由于初始He载气流速(100 mL/min)、进入分流接口时的流速(10 mL/min)与进入离子源时流速 (0.3 mL/min)的差异导致样品燃烧产生的目标气体中99.7%的气体被浪费, 样品的总体利用率仅有0.3%。因此如何减少样品在测试过程中的损耗, 提高样品的利用率, 从而减少需样量的关键在于缩小载气流速的差距。本实验在常规EA-IRMS测试技术基础上进行了关键改进,在元素分析仪和分流接口之间设计增加一个由六通阀和自动加热冷阱构成的装置, 自动加热冷阱可在Load模式时收集SO2气体, 六通阀在Load-Inject模式之间切换时可改变He载气流速, 通过与分流接口匹配的反吹He载气(10 mL/min)将冷阱中富集的SO2气体送入分流接口, 从而保证进入分流接口前样品燃烧产生的SO2气体全部收集。这一改进, 理论上可以将样品的利用率提高10倍, 该系统需硫量降至3~13 μg, 并且可以提高反应管的寿命, 降低清灰的频率, 提高工作效率。同时有效避免了拖尾的产生, 提高分析结果的精密度。本次微量样品实验获得的硫同位素数据与常规方法一致, 分析精度优于0.15‰(1SD), 测量值与真值的差异在0.4‰以内, 达到国际同类实验室先进水平。此外, 该方法可为微量有机碳、氮同位素EA-IRMS测试分析工作的开展提供参考经验。     

An Automatic Method for Detecting Transients and Variable Sources in AST3 Survey Based on Image Subtraction and Random Forest

AST3-2 (the second Antarctic Survey Telescope) is located in Antarctic Dome A, the loftiest ice dome on the Antarctic Plateau. It produces a huge amount of observational data which require a more efficient data reduction program to be developed. Also the data transmission in Antarctica is much difficult, thus it is necessary to perform data reduction and detect variable and transient sources remotely and automatically in Antarctica, but this attempt is restricted by the unsatisfactory performance of the low power consumption computer in Antarctica. For realizing this purpose, to develop a new method based on the existing image subtraction method and random forest algorithm, taking the AST3-2 2016 dataset as the test sample, becomes an alternative choice. This method performs image subtraction on the dataset, then applies the principle component analysis to extract the features of residual images. Random forest is used as a machine learning classifier, and in the test a recall rate of 97% is resulted for the positive sample. Our work has verified the feasibility and accuracy of this method, and finally found out a batch of candidates for variable stars in the AST3-2 2016 dataset.     

准噶尔盆地玛湖凹陷西斜坡三叠系百口泉组扇三角洲优质储集层预测

针对准噶尔盆地玛湖凹陷西斜坡风南地区三叠系百口泉组扇三角洲砂岩物性空间变化大、优质储集层(孔隙度大于7.4%,渗透率大于0.05×10^-3μm²)预测难的问题,在沉积岩石学、地震沉积学以及地震反演和解释理论指导下,综合利用测井、岩心和三维地震等资料开展了高精度层序地层划分、沉积微相描述和优质储集层地震反演研究。建立了风南井区四级层序地层格架,明确了扇三角洲多期水进水退的充填过程,指出SSQ3和SSQ5是优质储集层的发育层系;识别出扇三角洲平原分流河道、河道间和扇三角洲前缘水下分流水道、河口坝、席状砂等沉积微相,指出扇三角洲平原是优质储集层发育相带;通过应用高分辨层序地层纵向边界和沉积相横向边界约束,进行分层相控叠后地震波阻抗反演,提升储集层预测精度,在SSQ3和SSQ5预测5个优质储集层发育区,提出3口井的井位建议,钻探均获工业油流。     

深渊着陆器坐底位置的精确测量和反演计算

深渊着陆器在下潜过程中,由于受到海流的影响,其坐底位置会偏离投放点,偏离量和偏离方向很难估计,而声学水下定位方法在深渊海区不太适用.通过在海面进行气枪十字放炮作业,利用炮点准确位置和时间信息,以及海底地震仪(Ocean Bottom Seismometer,OBS)记录的直达水波到时,对深渊着陆器的坐底位置进行精确测量可以提供有效定位数据.该方法对采集得到的数据进行预处理,使用蒙托卡洛方法结合最小二乘法,对万泉号着陆器实际坐底位置进行了反演校正,得到三个潜次的位置校正偏移量分别为211m、178m、861m,偏移方向各不相同,校正精度为±20m,比以往的定位精度大为提高.此外,该方法测量确定了三个海底标志物的精确位置,为未来万米无人或载人潜器提供了可靠的水下参考坐标.     

陆地可控源电磁法探测效果的频率响应

陆地可控源电磁法的观测资料可依据频段范围近似地划分为近区场、中间区场及远区场,但采用测量相互正交电、磁分量,并计算视电阻率的资料处理方式只适用于远区场数据.为更有效地利用陆地可控源电磁法不同区间场的观测资料,本文结合三维数值模拟技术并采用电场分量直接进行反演的策略,对不同区间电场的响应特征与探测效果进行了分析.数值模拟结果表明:近区电场的异常响应最明显,异常响应不随频率发生显著变化,但纵向分辨能力差;远区电场异常响应随频率发生显著变化,其探测深度取决于频率的高低;中间区场较为复杂,地表电场异常响应的等值线中心并不是位于异常体中心正上方,而是在沿场传播方向上向异常体与围岩的分界面处偏移,并且发现中间区场资料的加入会影响反演结果的准确性.综合合成数据和野外实测资料的反演结果,发现结合近区场和远区场资料而舍弃中间区场资料的反演效果更佳,这为陆地可控源电磁法资料的反演解释提供了一种有效途径.     

Seismic amplitude inversion for the transversely isotropic media with vertical axis of symmetry

Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data.